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Abstract: 

The LLS stock market model is a model of heterogeneous quasi-rational investors 
operating in a complex environment about which they have incomplete information. We 
review the main features of this model and several of its extensions. We study the ef- 
fects of investor heterogeneity and show that predation, competition, or symbiosis may 
occur between different investor populations. The dynamics of the LLS model lead to 
the empirically observed Pareto wealth distribution. Many properties observed in actual 
markets appear as natural consequences of the LLS dynamics: truncated Levy distribution 
of short-term returns, excess volatility, a return autocorrelation "U-shape" pattern, and a 
positive correlation between volume and absolute returns. 

1. The LLS Model 

LLS is a microscopic representation model of the stock market. Its details and some 
generalizations of it can be found in [2] . In the present account we introduce the basic LLS 
ideas and the model main results. We consider below a market with only two investment 
options: a bond and a stock (see [3] for an extension to a multiple stocks case). The 
model involves a large number of virtual investors characterized each by a current wealth, 
portfolio structure, probability expectations and risk taking preferences. These personal 
characteristics come into play in each investor's decision making process as schematically 
seen in Fig. 1. 

The bond is assumed to be a risk-less asset yielding a return at the end of each time 
period. The bond is exogenous and investors can buy from it as much as they wish at a 
given rate. 

The stock is a risky asset with overall returns rate H[t) composed of two elements: 

(i) . Capital gain (loss): If an investor holds a stock, any rise (fall) in the market price 
of the stock contributes to an increase (decrease) in the investors' wealth. 

(ii) . Dividends: The company earns income and distributes dividends. 

Each investor i is confronted with a decision where the outcome is uncertain: which is 
the optimal fraction X{i) of his/her wealth to invest in stock? According to the standard 
theory of investment each investor is characterized by a utility function (of its wealth) 
U{W) that reflects his/her personal risk taking preference (here we take for simplicity 
U{W) = InW see [1] for a prospect theory extension). The optimal X{i) is the one that 
maximizes the expected value of his/her U{W) (we take into account all the unknown 
factors influencing decision-making (such as liquidity constraints or deviations from ratio- 
nality) by adding a small random variable (or "noise") to the optimal proportion X{i)). 
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The expected value of U (W) depends of course on the expected probabihties for the var- 
ious values of H to be realized in the future. In LLS the investors expectations for the 
future iif's are based on extrapolating the past values. More precisely each investor recalls 
the last k returns on the stock and expects that each of them may take place again with 
equal probability. The extrapolation range k differs between various investors and it will 
be in the sequel of this paper the main parameter inducing market inhomogeneity. 

At fixed time intervals each investor revises the composition of its portfolio and decides 
for a new market order. The aggregate of these orders determines the new stock price by 
the market clearance condition as explained below. Once each investor decides on the 
proportion of his/her wealth X{i) that (s)he wishes to hold in stocks, one can derive the 
number of stocks N{i,ph) it wishes to hold corresponding to each hypothetical stock price 
Ph- Since the total number of shares in the market N, is fixed there is a particular value 
of the price p for which the sum of the N{i,p) equals A^. This value p is the new market 
equilibrium price. Upon updating accordingly the traders' portfolios, wealth and list of 
last k returns, one is ready for the next market iteration. This process is repeated for each 
time step, and the market prices are recorded throughout the run. 

Figure 1: The Flow Chart of the LLS market framework 
2. Crashes, Booms and Cycles 

The LLS model provides already at the level of a quite homogenous traders population 
a convincing description of the emergence of cycles of booms and crashes in the stock 
markets. In a market with one species of investors all having a homogenous memory 
(extrapolation) range spanning the last k returns of the stock, the stock price alternates 
regularly between two very different price levels. The explanation for this behavior is as 
follows: 

Assume that the rate of return H{t) on the stock at a time t is higher than the oldest 
remembered return {H{t — k)). The addition of H[t) and the elimination of H{t — k) 
creates then a new distribution of past returns that is better than the previous one. Since 
the LLS extrapolating investors use the past k returns to estimate the distribution of the 
next period's return, they will be lead to be more optimistic and increase their investments 
in the stock. This, in turn, will cause the stock price to rise, which will generate an even 
higher return. 

This positive feedback loop stops only when investors reach the maximum investment 
proportion (i.e. X{i) = 100%: we do not allow borrowing or short selling), and can no 
longer increase their investment proportion in the stock. The dividend contribution to 
the returns is small compared with this high price at this stage. In the absence of noise 
the returns on the stock at this plateau converge to a constant growth rate which is just 
slightly higher than the riskless interest rate (see [4]). In other words, in the absence of 
noise the price remains almost constant, growing only because of the interest paid on the 
bond (more money entering the system and being invested in the stock). 

When there is some noise in the system the price fluctuates a little around the asymp- 
totic high level, because of the small random fluctuations in the investment proportions. 
These fluctuations generate some negative returns (on a downward fluctuation) and some 
high returns (when the price goes back up). One might suspect that a large downward 
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fluctuation might trigger a reverse positive feedback effect, trader expectations wiU lower, 
investment proportions will decrease, the price will drop, generating further negative re- 
turns and so on: a crash. This can happen during the "plateau" period but only after the 
previous sharp price boom which generated an extremely high return, is forgotten. And, 
indeed, this is exactly what happens. Since it takes k steps to forget the boom, the high 
price plateaus are a bit longer than the extrapolation span [k days to forget the boom 
+0(1) more days until a large enough negative fluctuation occurs). 

The crash generates a disastrous return and, until it is forgotten, investment propor- 
tions and hence the price remains very low. When the price is low, the dividend becomes 
signiflcant and the returns on the stock are relatively high (compared with the bond). 
Once the crash has been forgotten, all the returns that are remembered are therefore high, 
and the price jumps back up. Thus, the low price plateaus are k steps long. This completes 
one cycle, which is repeated throughout the run. This (quasi-)periodicity is best viewed 
in the Fourier transform of the price time evolution (Fig. 2) as a series of narrow peaks 
around the frequency 2k-\-0{l) and its harmonics (note however that the dynamics is not 
perfectly periodic and therefore in spite of its simplicity, according to some mathematical 
criteria it may fall into the "complex" category. In the present paper we reserve however 
the term "complex" for dynamics that are truly complicated to the degree that they do 
not admit simple verbal or mathematical description or understanding). 

The homogenous stock market described above exhibits booms and crashes. However, 
the homogeneity of investors leads to unrealistic periodicity. As shown below, when there 
is more than one investor species the dynamics becomes much more complex and realistic. 

Figure 2 : The Fourier transform of the price in a market with one species with 
extrapolation range A; = 10. The market contained 10000 traders that had initially equal 
wealth invested half in stock and half in bonds. 

3. Realistic Features in LLS with Many Species 

Our numerical experiments within the LLS framework have found that already a 
small number of trader species (characterized by different extrapolation ranges k) leads 
qualitatively to many of the empirically observed market phenomena. 

In reality, we would expect not just a few trader types, but rather an entire spectrum of 
investors. When the full spectrum of different trader species (fundamentalists and various 
other types - see [1] for the detailed operational deflnition) is considered it turns out that 
"more is different" [5]: the price dynamics becomes realistic: booms and crashes are not 
periodic or predictable, and they are also less frequent and dramatic. At the same time, 
we still obtain many of the usual market anomalies described by the experimental studies 
(however in the limit of inflnite times or inflnite number of investors, the dynamics may 
revert to predictable patterns [3]). 

We list below a few such realistic features: 

Return Autocorrelations: Momentum and Mean-Reversion 

In the heterogeneous population LLS model trends are generated by the same positive 
feedback mechanism that generated cycles in the homogeneous case (section 2): high (low) 
returns tend to make the extrapolating investors more (less) aggressive, this generates 
more high (low) returns, etc. 
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The difference between the two cases is that in the heterogeneous case there is a very 
complicated interaction between all the different investor species and as a result there are 
no distinct regular cycles but rather, smoother and more irregular trends. There is no single 
cycle length - the dynamics is a combination of many different cycles corresponding to the 
many extrapolation ranges k. This makes the autocorrelation pattern also smoother and 
more continuous. The return autocorrelations in the heterogeneous LLS model conform 
to the empirical findings: In the short run the autocorrelation is positive - this is the 
empirically documented phenomenon known as momentum: high returns during a trading 
quarter tend to be followed by more high returns in the following months, (and low returns 
tend to be followed by more low returns). In the longer run the autocorrelation is negative 
(after a few years of boom, one usually experiences a few "dry" years), which is known as 
mean-reversion. For even longer lags the autocorrelation eventually tends to zero [1]. The 
short run momentum, longer run mean-reversion, and eventual diminishing autocorrelation 
creates the general "U-shape" that is found in empirical studies [7]. 

Excess Volatility 

In markets with a large fundamentalist population (see [1] for their detailed operative 
definition in the LLS model), the price level is generally determined by the fundamen- 
tal value of the stock. However, the market extrapolating investors occasionally induce 
temporary departures of the price away from the fundamental value. These temporary 
departures from the fundamental value make the price more volatile than the fundamental 
value. 

Following Shillcr's [8] methodology we measured the standard deviations of the de- 
trended price and fundamental value. Averaging over 100 independent simulations we 
found [1] respectively 27.1 and 19.2, which is an excess volatility of 41%. 

Heavy Trading Volume 

In an LLS market with both fundamentalists and market extrapolating investors (over 
various k ranges), shares change hands continuously between the various groups: 

When a "boom" starts, the extrapolating investors observe higher ex-post returns and 
become more optimistic, while the fundamentalists view the stock as becoming overpriced 
and become more pessimistic. Thus, at this stage the market extrapolators buy most of 
the shares from the fundamentalists. 

When the stock crashes, the opposite is true: the extrapolators are very pessimistic, 
but the fundamentalists buy the stock once it falls below the fundamental value. Thus, 
there is substantial trading volume in this market. The average trading volume in a typical 
LLS simulation was about 1, 000 shares per period, or about 10% of the total outstanding 
shares. 

Volume is Positively Correlated with Absolute Returns 

The typical scenario in an LLS run is that when a positive trend is induced by the 
extrapolating investors, the opinions of the fundamentalists and the extrapolating investors 
change in opposite directions: 

- The extrapolating investors see a trend of rising prices as a positive indication about 
the future return distribution, while 

- The fundamentalists believe that the higher the price level is (the more overpriced 
the stock is), the harder it will eventually fall. 
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The exact opposite holds for a trend of faUing prices. Thus, price trends are typicaUy 
interpreted differently by the two investor types, and therefore induce heavy trading vol- 
ume. The more pronounced the trend (large price changes), the more likely it is to lead 
to heavy volume. 

In order to verify this relationship quantitatively we regressed volume V{t) on the 
absolute rates of return r{t) for 100 independent simulations. We run the regressions: 

V{t) — a + b\r{t) \ + random{t) 

We found an average value of 870 for h with an average t- value of 5.0. Similar results were 
obtained for time lagged-returns. 

4. Predation, Competition and Symbiosis between Trader Species 

In section 2 it was explained that a homogenous population of traders that extrapolate 
the last k returns leads to cycles of booms and crashes of period 2A;-|-0(1). When there are 
two species with extrapolation ranges ki and respectively /c2, we observe sharp irregular 
transitions between eras where one species dominates (cycles of period 2ki + 0(1)) and 
market eras where the other species dominates (cycles of period 2k2 + 0{1)). When the 
number of trader species is three, there are dramatic qualitative changes: generically, the 
dynamics becomes complex. We show that complexity is an intrinsic property of the stock 
market. This suggests an alternative explanation to the widely accepted but empirically 
questionable random walk hypothesis. We discuss below some of the market ecologies 
possible with only two species of traders. Of course the picture becomes more complex 
later, when 3 or more species are introduced. 

Market Ecologies with Two Trader Species 

When there are two trader species with different extrapolation spans it turns out that 
the nature of the dynamics is determined by the ratio of the extrapolation spans of the 
two species. In [9] we performed a qualitative theoretical analysis of this phenomenon 
and supported it by microscopic simulations. We showed that in market eras in which 
one species (of extrapolation range ko) dictates the dynamics (i.e. boom-crash cycles have 
periods of length 2/co+0(l)) the second species (with extrapolation range k) has generically 
the following performance: 

A : If /cq < < 2/co then k is performing very poorly (looses money) 
B : If 2nfco < k < (2n + l)fco (with n natural number), then k is doing relatively well 
C : If (2?i + l)/co < k < 2n/co, n > 1, then k does better than in A but worse than in B 
D : If /c < /cq, then k is doing well . 

These facts turned out sufficient to understand the main 3 cases that a 2-species 
ecology can display: 

Case 1: predator - prey dynamics 

If one considers one species with an extrapolation range ki = 10 and a second species 
with an extrapolation range /c2 = 14 it turns out that the resulting ecology dynamics is a 
predator-prey one. In fact the LLS market dynamics leads in this case to the extinction 
(total impoverishment) of the ki species: after some time the entire wealth on the market 
belongs to the species ki = 10 (Fig. 3). As a consequence, the market price presents clear 
cycles of booms and crashes of periodicity clustered around 24 = 2*10-1-0(1). 
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This is easily understood since according to the property A above the ^2 = 14 popula- 
tion is performing poorly when the ki dictates the market periodicity while the population 
10 is performing well according D in the hypothetical periods when k2 = 14 dictates the 
market periodicity. This is only an example of a large class of parameters that lead to 
predator- prey systems and which may result in the total extinction of one of the species. 

Figure 3 : Fraction of the wealth that the species ki = 10 possesses in Case 1. The 
traders in the market belonged to 2 species consisting each of 5000 traders. Each trader 
owed at the beginning 5000 dollars in cache and 5000 shares (worth each I.4 dolars). 

Case 2: competitive species 

If one chooses ki = 10, k2 = 26, the species with extrapolation range 26 gains during 
the periods when the species ki = 10 dominates (property B) but species ki = 10 gains 
when the species k2 = 26 dominates (property D). It is therefore reasonable that one 
species can not dominate the other indefinitely. Indeed, a look at the fraction of the 
wealth held by the species with extrapolation range A;i = 10 reveals alternating eras of 
dominance (Figure 4). This is also reflected in the alternance between price cycles (~ 56) 
corresponding to /c2 = 26 and price cycles (~ 24) corresponding to ki = 10. Clearly this 
alternance between the 2 species corresponds to a classical competitive ecology, in which 
two competing species take turns in dominating the ecology. Note however that most of 
the time it is the population k2 which dominates the wealth. This seems to be a generic 
tendency in the long runs limit. 

Figure 4- Fraction of the wealth that the species ki = 10 possesses in Case 2. The 
initial conditions were similar to Figure 3. 

Case 3 symbiotic species 

In the case ki = 10, /c2 = 36, similarly to the 10 — 26 market, the investors with 
extrapolation range /c2 = 36 are doing better than those with extrapolation range ki = 10 
when ki — 10 dictates the dynamics (cf. property C). On the other hand ki = 10 are doing 
better when the species /c2 = 36 dictates the dynamics (cf. D). Hence, we may speculate 
that again we will find alternating eras of dominance. Figure 5 shows that this is not the 
case. The difference between this case and the 10 — 26 case is that here the market remains 
stuck in a "metastable" state: the extrapolation range 36 population never gains enough 
wealth to dictate long cycles. Thus, the system remains in a state of symbiosis throughout 
the run: the price cycles correspond to the short species extrapolation range span ki — 10 
while 70 — 80% of the wealth stays with the long extrapolation span species k2. 

For very long k2 extrapolation ranges, the share of the total wealth detained by k2 
can be even larger (approaching unity). 

Figure 5 -Fraction of the wealth that the species ki — 10 possesses in case 3. 

In conclusion [9] has uncovered a quite lively ecology of the traders populations in the 
LLS model and 
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"observed phenomena ranging from complete dominance of one population to alternat- 
ing eras of domination and to symbiosis. . . . Our results suggest that complexity is an 
intrinsic property of the stock market. The dynamic and complex behavior of the market 
need not be explained as an effect of external random information. It is a natural prop- 
erty of the market, emerging from the strong nonlinear interaction between the different 
investor subgroups of the market ..." 

The main source of endogenous dynamics in the LLS model turns out to be the 
feedback between the market price fluctuations and the wealth of the investors belonging 
to various species: 

- On one hand the wealth of the investors determines their influence on the price 
changes (at the short range): e.g. the richest determine the periodicity of the boom-crash 
cycles. 

- On the other hand, the variations in the price determine changes in the distribution 
of wealth, which iterated over longer time intervals, result in changes in the market price 
cycle periodicity regime. 

The entire cycle of rise and fall of a given species can be schematically described as: — > 
The species has by chance a (momentary) winning strategy — > Investors belonging to the 
species gain wealth Overall wealth of the investors belonging to the species increases — > 
Bids of investors belonging to the species become large Investor bids influence the market 
price adversely (self-defeating) — > Trading of investors belonging to the species becomes 
ineflBcient — > Investors lose money — > Investors belonging to the species become poor 
— > Species wealth and market relevance vanish Other species with different strategies 
become winners Cycle re-starts (with the new winning strategies). 

A few comments are in order: 

1. the concept of efficient strategy is only a temporary one as it depends crucially of 
the state of the market: by its very efliciency at a certain moment, a strategy prepares the 
seeds of its failure in the future. 

2. the biological and cognitive analogies are useful but their limits should be under- 
stood: 

- in biology, the species selection mechanism is based on the disappearance of the 
inefficient individuals. 

- in the learning adaptive agents' case, the individuals discard loosing strategies for 
new ones. 

In the LLS market framework, while it is possible to include the above effects, they are 
not necessary: the strategies selection takes place automatically by their carriers (traders 
belonging to the species) losing or gaining: for the market to be efficient, no a priori 
intelligence nor explicit criteria for the evaluation and comparison of market performance 
are required: just the natural (Adam Smith's "invisible hand") market mechanisms. 

3. While the adverse influence on the market price implied by the large orders coming 
from rich agents' leads automatically to inefficiency in their operations (except for rich 
agents which follow a buy-and-hold strategy and therefore do not influence (adversely) the 
market), the mere lack of market influence due to poverty does not guarantee a winning 
strategy. It is necessary therefore that there are enough strategies and enough agents in 
the market for insuring its efficiency. 
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Three Investor Species 

One might suspect that the three species dynamics is a natural extension of the two 
species dynamics. Instead of alternating between two cycle lengths the system may just 
alternate between the three possible states of dominance. Figures 6-7 shows that this is 
not at all the case. These figures depict a typical part of the dynamics of a three species 
market, with extrapolation spans of 10, 141 and 256 respectively. With the introduction 
of a third species the system has underwent a qualitative change: there is no specific 
cycle length describing the time series. Instead, we see a mixture of different time scales: 
the system has become complex. Prediction becomes very difficult, and in this sense the 
market is much more realistic. Figure 6 shows the power struggle between the three species 
while the Figure 7 depicts the Fourier transform of the price evolution during this run. 

Although the dynamics is complex, it is clear from Figures 6 and 7 that there is an 
underlying structure, which perhaps may be analyzed by the properties A, B, C, D and 
their generalizations. For instance it would appear from the Figure 6 that 141 and 256 take 
turns in dominating while 10 has a chance to a non- vanishing wealth share only occasionally 
in the transition intervals between 141 and 256 dominated eras. The dynamics generated 
by only three investor species can be extremely complex, even without any external random 
influences. 

Figure 6: The species wealths in a market with 3 species of extrapolation ranges of 

respectively 10, I4I and 256 days. Initially the 3 species possessed equal wealth distributed 
equally between stock and bond. Each species consisted of 1000 traders. 

Figure 1: Fourier transform of the stock price time evolution in the market described 
in Fig 6. 

5. Generalized Lotka-Volterra models for markets with multiple species 

Inspired by the above facts we devised an effective dynamics that stylized the features 
uncovered in the LLS model and extended them to a more generic framework. Instead 
of following in detail the way the market price influence each species and individual i , 
we assumed that this influence can be represented through multiplying their wealth Wi{t) 
by stochastic multiplicative factors \i{t). This is natural in the LLS model in which 
the investments of the individuals (and consequently their returns) are fractions of their 
wealth (as implied by the constant relative risk aversion utility functions). The stochastic 
proportionality between personal returns and personal wealth is consistent with the real 
data that show that the (annual) individual income distribution is proportional to the 
individual wealth distribution [10]. 

We proposed [11] therefore a model including the above stochastic autocatalytic prop- 
erties of the capital as well as the cooperative, diffusive and competitive/ predatory in- 
teractions between the species. The resulting model turned out to be a straightforward 
generalization [11] of the Lotka-Volterra system (discrete logistic equation) well known 
previously in population biology: 

Wi{t + 1) = Xi{t)wi{t) + '^akWk{t) - Wi^bkWkit) 

k k 
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Where the sum is over all the N traders participating in the market. There are a few 
other (mutually non-exclusive) possible interpretations for Wi in addition to the individual 
wealth: the wealth associated with a particular investing strategy, the capitalization asso- 
ciated with a particular company /industry, or the number of investors following a common 
trend (herd). 

Particularly interesting cases were studied subsequently: 

1) The linear case where the total wealth diverges to larger and larger values (inflation, 
production) : 

Wi{t + 1) = \i{t)wi{t) + ^ akWk{t) 

k 

2) The case in which the binary interactions between individuals are expressible in 
terms of interactions with the total wealth W = J2k "^k'- 

Wiit + 1) = Xi(t)wi{t) + aWit) - bwiWit) 

3) The case in which individual wealth is bounded from below by a certain fraction c 
of the average wealth w = W/N [12]: 

Wi{t + 1) = Xi{t)wi{t) 

except if 

Xi{t)wi{t) < cw 

when 

Wi{t + 1) = cw 

4) The case of the random multiplicative wealth dynamics 

Wi{t + 1) ^ Xi{t)wi{t) 

with variable number of traders [13]: 

- traders which fall below a certain fixed minimal wealth Wmin drop from the market 

- a number of traders proportional to the total wealth increase: 

AN = ciW{t + 1) - Wit))/Wmin 

join the market at each time step i.e. the average wealth remains constant in this process: 

w = Wmin/c 

One assumes that each of the new traders brings an initial investment equal to Wmin 
which means that the total amount of added wealth is 

ANWmin=c{W{t+l)-W{t)) 

i.e. a fraction c of the total wealth increase {W{t + 1) — W{t)). 
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In most of these systems were assumed asynchronous: at each time step, only one 
(randomly chosen) Wi was updated. Very striking generic results can be obtained with all 
these models in certain relevant regimes. We limit ourselves below to the universal scaling 
properties (power laws). 

6. Pareto Law in LLS and Lotka-Volterra models 

The efficient market hypothesis and the Pareto law are some of the most striking and 
basic concepts in economic thinking. It is therefore very important that our models above 
succeed to connect them in a very essential way. 

Let us discuss this in more detail. More than a hundred years ago, Pareto [14] dis- 
covered that the number of individuals with wealth (or incomes) with a certain value w 
is proportional to This later became known as the Pareto Law. The LLS model 

treats the individual investor wealth as a crucial quantity, and it views its feedback rela- 
tion to the market dynamics as the main source driving the endogenous dynamics of the 
market. 

It turns out that in the conditions in which the participants in the market do not 
have a systematic advantage one over the other (which is in fact expected in an efficient 
market), a dynamics of the LLS type leads always to a Pareto law. The actual value of the 
exponent a depends on the particular parameters used in the model. Mainly, as explained 
below a is influenced by the social security policy. If one does do not allow any individual 
to become poorer than a certain fraction c of the current average wealth then, for a wide 
range of conditions a = 1/ ( 1 — c) . This is confirmed in Figure 8 which plots the wealth 
distribution in the LLS model with k — 3 and c = 0.2 (and U{W) — InW). 

Figure 8: The wealth distribution of the investors in an LLS model with a poverty line 
of c = 20% of the average wealth. One a double logarithmic scale one obtains a straight 
line with slope 2.2 corresponding to an a of 1.2. 

The market consisted of 10000 traders and the measurement was performed as a "snap- 
shot" after 1 000 000 "thermalization" market steps. Initially all the traders had equal 
wealth ($1000) equally distributed between bond and stock. 

In fact it has been proven [11-13] theoretically that any of the effective dynamics of 
the type 1-4 with Aj distribution independent on i or Wi leads always to a power law Pareto 
distribution. 

In a wide range of models, the generic rule is that 
a = l/(l-c) 

where c is essentially the market global impact factor [13]: 

c= exogenous new capital ADDED to the market / increase in stock cap- 
italization due to market price increase 

Since the increase in the capitalization is the increase in wealth that the owners incurr 
upon their investment of new capital, the ratio can be also interpreted as the long range 
market return factor 

c = 1/ (long range market return factor). 

Let us explain in short how such results were obtained [15]. The crucial observation 
is that for large Wi values, the non-stationary multiplicative system of interacting tUj's is 
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formally equivalent to a statistical mechanical (additive) system in thermal equilibrium 
when expressed in terms of the variables Ui{t) = lnwi/w{t). For instance the system 3) is 
mapped into a system of particles diffusing in an energy potential field u with a ground 
level uq — In c. In thermal equilibrium, all such systems (independent on the details of the 
interactions between their particles) have an universal probability distribution discovered 
by Boltzmann more than 100 years ago: 

P{u) ~ exp{—au) 

When re-expressed in terms of the original Wi variables this gives a Pareto power law 
distribution: 

P{w) ~ w-^-'^ 

The exponent a can be estimated from the integrals representing the total wealth and the 
total number of traders [12]. For instance, in the models 3-4, in the limit of N ^ oo , the 
result is: 

a = 1/(1 - Wmin/w) 

i.e. 

a = 1/(1 -c) 

Thus, the Pareto law is the exact analogue of the Boltzmann law for stochastic systems 
that are multiplicative rather then additive. 

For finite N the a is given by the implicit transcendental equation [16]: 

N = [((1 - {N/cr)/a]/[{{l - {N/cr-')/{a - 1)] 
which for N <^ e^/^ gives approximately: 

a - lniV/(ln(iV/c)) < 1 

which incidentally means that in this regime all the wealth belongs to only a few individuals. 

In the system 1 defined above, in the appropriate thermodynamic limit The analog 
result is: 

a= 1/(1 -c) 

with 

c~ 2a/(< A > +cr^/2) 

where a is the standard deviation of A 
And in the case 2: 

cr^2a/{a^ + a'^) 

independent on b and < A >. 

7. Market Efficiency, Pareto Law and Thermal Equilibrium 

The formal equivalence between the non-stationary systems of interacting tUj's and the 
equilibrium statistical mechanics systems governed by the universal Boltzmann distribution 
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has far reaching imphcations: it relates the Pareto distribution to the efficient market 
hypothesis: In order to obtain a Pareto power law wealth distribution it is necessary and 
sufficient that the returns of all the strategies practiced in the market are stochastically 
the same, i.e. there are no investors that can obtain "abnormal" returns. 

Therefore, the presence of a Pareto wealth distribution is a measure of the market 
efficiency in analogy to the Boltzmann distribution whose presence is a measure to thermal 
equilibrium. Indeed physical systems which are not in thermal equilibrium (e.g. are forced 
by some external field - say by laser pumping) do not fulfill the Boltzmann law. Similarly, 
markets that are not efficient (e.g. when some groups of investors make systematically 
more profit than others) do not yield power laws (see Fig 9). Optimal market and power 
laws are the short time and long time faces of the same medal/phenomenon. 

This analogy is consistent with the interpretation of market efficiency as analog to the 
Second law of Thermodynamics: 

- one can extract energy (only) from systems that are not in thermal equilibrium 

- one can extract wealth (only) from markets that are not efficient. 

- by extracting energy from a non-equilibrium thermal system one gets it closer to an 
equilibrium one. 

- by extracting wealth from a non-efficient market one brings it closer to an efficient 

one 

-in the process of approaching thermal equilibrium, one also approaches the Boltzmann 
energy distribution 

- in the process of approaching the efficient market one also approaches the Pareto 
wealth distribution. 

-by having additional knowledge on a thermodynamic system state one can extract 
additional energy (e.g. Maxwell demons gedanken experiment) 

-by having additional knowledge on a financial system one can extract additional 
wealth. 

This double analogy 

thermodynamic equilibrium ~ efficient market 
Boltzmann law ~ Pareto law 
holds in the details of their microscopic origins: 

- the convergence to statistical mechanics equilibrium depends on the balance of the 
probability fiow entering and exiting each energy level. This is usually insured microscop- 
ically by the fact that the a priori probability for a molecule to gain or loose an energy 
quanta in a collision is the same for any energy level with the exception of the collisions 
including molecules in the ground state which can only receive (but not give) energy. 

- in the stochastic models 1-4, the convergence of the wealth to the power-law is 
insured by the balance of fiow of investors from one level of [log (relative wealth)] to 
another. At the individual level, this is enforced by all the individuals having the same 
(relative) returns probability distribution (except for the individuals possessing the lowest 
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allowed wealth). If this condition is not fulfilled, one does not get a wealth distribution 
power law. 

These facts should guide us in the practical runs in establishing which combinations 
of strategies (or the strategy selection strategies) are producing a realistic market "in the 
Pareto sense". In Figure 9 one sees the wealth distribution in a model in which there 
are 2 trader species with slightly different distributions of A. One sees that even a small 
violation of the A uniformity leads to significant departures from the Pareto law which are 
inconsistent with the historical experimental facts. The absence of such departures in real 
life is a strong indication of the market efficiency in the weak stochastic sense (that all 
investors have stochastically the same relative returns distribution). 

Figure 9: Wealth distribution for 2 investor species with different return distributions. 
Model 3 was used with a lower wealth bound of c — 20%. 

A is randomly drawn. For the first species X is 1.10 or 0.95 with equal probability. For 
the second "more talented" species A is 1.11 or 0.96 with equal probability. The 2 species 
were each composed of 10000 traders with initially equal wealth (1000 dollars each). 

The measurement of the wealth distribution was performed after a "thermalization 
period" of 100 000 wealth updatings. 

8. Leptokurtic Market Returns in LLS 

It has been long known that the distribution of stock returns is leptokurtic or "fat- 
tailed". Furthermore, a specific functional form has been suggested for the short-term 
return distribution (at least in a certain finite range) - the Levy distribution [17]. This 
feature is present in the LLS model, and is directly related to the Pareto distribution of 
wealth. 

The central limit theorem insures that in a wide range of conditions the distance 
reached by a random walk of t steps of average squared size is a Gaussian with standard 
deviation s\/t. Suppose that at time t — Q one has N positive numbers Wi(0); i = 1, A?" 
of order 1 and sum 1^(0). Suppose that at each time step one of the numbers varies 
(increases or decreases) by a fraction s j {t) <^ 1 extracted from a random distribution with 
average squared (and mean). What will be the probability distribution of the sum 
W{t) after t steps? According to the central limit theorem this would be the Gaussian 

P{W,t) = l/(V2^)e-(^W-^(°))'/2*^' 

since it consists of t steps of average squared size s^. 

If one interprets Wi{t) as the value of the stocks owned by the trader i at time t, then 
W{t) — '^^Wi is the total market value of the stock and therefore {W{t) — W{0))/W{0) 
is the relative stock return for the time interval t. 

One sees that if the central limit theorem would hold, one would predict a Gaussian 
stock returns distribution. This is in fact the case for real stocks and time intervals longer 
than a few weeks. For significantly shorter times t however, the distribution of returns is 
very different from a Gaussian. Even though the exact shape of the returns distribution is 
not yet established experimentally, it is generally agreed that in certain ranges (typically 
"in the tails"- i.e. for large Wi values) it fits better a power law rather than a Gaussian. 
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Such a situation can in principle be explained by the following scenario: 
Suppose that at time t = one has an arbitrarily large number of positive numbers 
Wi{0). Suppose moreover that the probability distribution for the sizes of Wj(0) is 

P{w) ~ -u;-^-" 

Suppose that at each time step one of the Wj's varies (increases or decreases) by a fraction 
Si{t) <^ 1 of average squared size s^. What will be the probability distribution of the 
variation of the tu^'s sum W{t) — W{0) after t steps? 

One is tempted to think that the correct answer is given by 

P{W,t) = l/(y2^)e-(^W-^(°))'/2*^' 

for some s. However this is wrong. Indeed, assuming such an s exists would imply that 
the probability for the sum variation W{t) — W{0) to be 10 after a time t = l/(2s^) is: 

P{W{t) = W{0) + 10,t = l/{2s^)) ~ e-^°' ~ 10-32 

while in reality a lower bound for the probability of getting W{t) — W{0) = 10 in just one 
step it is obviously that given by 

P(w) ~ -u;-^-" 

I.e. P{W{t) = W{0) + 10, t = 1/(25^)) is at least of order 10"^"" which for a < 2 means 
it is larger than 10"^ ! 

This coarse estimations highlights the difference between the Gaussian distributions 
and the distributions generated by random walks with power distributed step sizes (called 
Levy distributions [18,17]): the presence of Wi^s of arbitrary size implied by a power law 
distribution insures that the large returns distribution is dominated by the power law of 
the individual step sizes rather than the combinatorics of the multiple events characterizing 
the Gaussian system. 

One sees now that the systems 1-4 (and consequently LLS) are exactly of the type 
one needs to explain returns distribution power tails: 

- on one hand according to section 6, the models 1-4 (and consequently LLS) insure a 
power distribution of Wj's. 

- on the other hand, in the models 1-4 the variation of the stock index W(t) is the 
sum of the variations of the individual Wi{t)^s. 

- these variations Wi{t + 1) — Wiit) are stochastic fractions siit) = Xi{t) — 1 of Wi 
as above (the fact that Xi{t) — 1 has not mean is taken care by working actually with 
Ui = \n{wi/id)). 

Therefore, according to the argument above, the effective models 1-4, which reflect 
the stochastic proportionality in LLS between individual wealth, individual investments 
and individual gains/losses predict that the price fluctuations in the LLS model will obey 
a Levy distribution (and in particular flt a power in some range of the "tail"). 

There is a proviso for this argument to hold: the number of individual terms N has 
to be larger than the number of time steps t. Otherwise the flnite size of the sample of 
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tUj's will show up in the absence of sizes Wi larger than a certain value. In fact for t N one 
recovers (slowly) the Gaussian distribution. 

In the LLS case, if the portfolio updatings are performed simultaneously by all the 
investors, the unit time step corresponds already to a time t — N. In order to verify the 
(truncated) Levy distribution and the power "tail" predictions, one has to look at the 
dynamics at a finer time scale. We therefore performed [12] LLS runs in which at each 
time step only one trader i reconsiders its portfolio investment proportion X{i). In such 
conditions, one expects to obtain a distribution which fits in a significant range a power 
law (up to large Wi values where the finite effects become important). 

This is in fact confirmed by the numerical experiments. While for the global updating 
steps one gets a Gaussian distribution, for the trader-by-trader procedure one obtains a 
truncated Levy distribution (Fig. 10). 

Figure 10: The returns distribution in the LLS model in which only one trader re- 
evaluates his/her portfolio per unit time, c = 0.2, k = 2t, U= \nW. 

The market contained 10 000 traders with initially equal wealth and portfolio compo- 
sition (half in stock and half in bonds). 

The number of market returns in intervals of 0. 001 were measured during 5 000 000 
market steps (after an initial 1 000 000 equilibration period). 

Note that in the central region of the short time returns (before the cut-off becomes 
relevant) the Levy distribution is characterized by an a equal to the exponent a of the 
traders' wealth distribution. 

As explained in Section 6, in certain conditions (e.g. model 4) one can interpret a as 

a = 1/(1- l/(long term market return factor)). 

Therefore the analysis above relates the stochastic distribution of the short term 
returns to the value of the long term returns via the exponent of the Pareto power 
law of individual incomes /wealths. 

Moreover the long term returns are related (e.g. model 3) via the value of the 

Pareto exponent a to the ratio (w/wmin) between the average wealth/income and the 
lowest admissible wealth/income: the value cu ~ 1.4 implies (cf. models 3-4) for both 
these quantities values of the order of 

1/c = a/{a — 1) 3.5. 

Speculatively, one may try to use the above relation in order to explain the stability of 
the Pareto constant a over the past century (and over the various countries and economies). 

Indeed one may relate the implied value 3.5 for both w/wmin and the long term market 
return to some basic biological invariant which is the average number of dependents / 
offsprings humans have: 

- if Wmin is the minimal amount necessary to keep alive one person in a certain society 
(cost of life), then the average income w will have to equal roughly Wmin times the number 
of dependents the average household head has to support. 

- at the social level, the total effort/wealth that one generation invests in the economy 
has to result in an economical growth capable to support a population larger by a factor 
equal to the average number of descendents. 
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Figure Captions 

Figure 1: The Flow Chart of the LLS market framework 

Figure 2 : The Fourier transform of the price in a market with one species with extrap- 
olation range k = 10. The market contained 10000 traders that had initially equal wealth 
invested half in stock and half in bonds. 

Figure 3 : Fraction of the wealth that the species ki = 10 possesses in Case 1. The traders 
in the market belonged to 2 species consisting each of 5000 traders. Each trader owed at 
the beginning 5000 dollars in cache and 5000 shares (worth each I.4 dolars). 

Figure 4- Fraction of the wealth that the species ki = 10 possesses in Case 2. The initial 
conditions were similar to Figure 3. 

Figure 5 -Fraction of the wealth that the species ki = 10 possesses in case 3. 



Figure 6: The species wealths in a market with 3 species of extrapolation ranges of re- 
spectively 10, 141 and 256 days. Initially the 3 species possessed equal wealth distributed 
equally between stock and bond. Each species consisted of 1000 traders. 
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Figure 7: Fourier transform of the stock price time evolution in the market described in 
Fig 6. 

Figure 8: The wealth distribution of the investors in an LLS model with a poverty line of 
c = 20% of the average wealth. One a double logarithmic scale one obtains a straight line 
with slope 2.2 corresponding to an a of 1.2. 

The market consisted of 10000 traders and the measurement was performed as a "snapshot" 
after 1 000 000 "thermalization" market steps. Initially all the traders had equal wealth 
($1000) equally distributed between bond and stock. 

Figure 9: Wealth distribution for 2 investor species with different return distributions. 

Model 3 was used with a lower wealth hound of c = 20%. 

A is randomly drawn. For the first species A is 1.10 or 0.95 with equal probability. For the 
second "more talented" species \ is 1.11 or 0.96 with equal probability. The 2 species were 
each composed of 10000 traders with initially equal wealth (1000 dollars each). 
The measurement of the wealth distribution was performed after a "thermalization period" 
of 100 000 wealth updatings. 

Figure 10: The returns distribution in the LLS model in which only one trader re-evaluates 
his/her portfolio per unit time, c = 0.2, k = 3, U= \nW. 

The market contained 10 000 traders with initially equal wealth and portfolio composition 

(half in stock and half in bonds). 

The number of market returns in intervals of 0. 001 were measured during 5 000 000 market 
steps (after an initial 1 000 000 equilibration period). 
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